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Abstract. We give an elementary construction of a time-invertible Markov 
process which is discrete except at one instance. The process is one of the 
quadratic harnesses studied in and It can be regarded as a random 

joint of two independent Poisson processes. 

1. Introduction 

According to '11", a stochastic process (Xf)j>o has the time inversion property, 
if it has the same finite-dimensional distributions as the process (iXx/f)t>o- In 
papers and the authors give criteria for the time-invertibiHty of Markov 
processes with transition probabihties that have smooth densities with respect to 
the Lebesgue measure. 

In this note we give a new elementary example of a time-invertible Markov 
process for which all transitions except to time t = 1 are discrete, see Proposition 
14.21 This improves upon ^ Corollary 3.4], where we gave a less elementary example 
of a time-invertible Markov process that had transition probabilities with a discrete 
component. Both examples are particular cases of a more general family of Markov 
processes which in 2 we called the bi-Poisson processes. According to 2, Example 
4.8 and Proposition 4.13], a bi-Poisson process with parameters (rj, 9, q) is a square- 
integrable Markov process {Xt ) which is uniquely determined by the following three 
properties: 

(1) E(Xt) = 0, E{XtXs) = mm{t, s}, 



11 — f f — S 

(2) E{Xt\Ts,u) = X, + - 

u ~ s u ~ s 

(3) Var(Xt|.F,,„) = 

{u -t){t- s) ( uX, - sXu ^ Xu-X, {uXs - sX„) (X„ - X, 

1 + rj 1- (1 - q)- 



u — qs \ u ~ s u — s {u — s) 

for all < s < < < u, where 

^s,u — o'iXt : < t < s or t > u}. 
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Property (O says that the bi-Poisson process is a harness, see [O]- Condition Q 
means that it is a quadratic harnesses, see [2]. The adjective "classical" refers 
to the value of parameter q = 1, compare |21 Section 4.2] and 'F, Section 4.2]. 
In Proposition 14. II we show that (Xt) can be constructed by joining together two 
independent Poisson processes with the same random gamma intensity. This is 
accomplished by appropriate affine transformations and deterministic changes of 
time. 

In 0] we use orthogonal polynomials to construct the transition probabilities of 
the bi-Poisson process when q = 0, and we show that its univariate distributions 
form a semigroup with respect to a certain generalized convolution related to free 
probability. The univariate distributions of the general bi-Poisson process were 
implicitly identified in [21 Example 4.8], and the corresponding Markov process is 
under construction in . This construction relies heavily on cumbersome identities 
between certain multi-parameter families of orthogonal polynomials, and identifies 
the transition probabilities in implicit form only. However, when q = 1, the explicit 
transitions probabilities can be read out. They turn out to be related to the pure 
birth and the pure death processes and are amenable to explicit elementary analysis. 

In this paper we present an elementary construction of the bi-Poisson process 
with parameters {r],9,l). Throughout most of the paper the value of the third 
parameter is fixed as g = 1, in which case we say that we consider a bi-Poisson 
process with parameters (rj, 9), skipping the third parameter of the triple. 

The plan of the paper is as follows. In Section |21 we give the construction of the 
process. In Section|21we verify that the construction indeed gives a bi-Poisson pro- 
cess. In Section0]we deduce some additional properties, including time-invcrtibility. 

2. Construction 

It is known, see |21 Proposition 4.13] that a bi-Poisson process (Xt) with param- 
eters (?7, 9) satisfies r]9 > 0. 

In the degenerate case 6*77 = it is known that the bi-Poisson process (Xt) 
is cither Xt = Bt, where (Bt) is the standard Brownian motion for rj = 9 = 0, 
or Xt ~ 9Ntfg2 — t/9, where (Nt) is the standard Poisson process when 77 = 0, 
9 Q, see J3 Theorem 1]. Passing to the time inverse (tXi^t), we see that 
Xt — •qtNi/(j^^2-^ — l/rj in the remaining degenerate case = 0, ry 7^ 0. 

We will therefore concentrate on the case 9r] > 0. Passing to {-—Xt) preserves (PJ 
replacing parameters (77, 9) by {—rj, —9), so we may assume 7;, 6* > 0. Replacing pro- 
cess (Xt) by process {y/7]j9Xtg/j,), we get the bi-Poisson process with parameters 
(v^, y/r]9). Thus without loss of generality we may assume that rj = 9 > 0. 

The moment of time i = 1 is preserved by the time-inversion and plays a special 
role in the construction. The bi-Poisson process traverses a family of deterministic 
lines, with jumps in the upwards direction when t < 1 and in the downwards 
direction when t > 1, see Fig. ^ 

The process is determined by specifying an integer that describes the line being 
followed at time t. The integers that describe the upwards jumps form a linear pure 
birth process with immigration with the time transformed to run on the interval 
[0,1). At time t = 1 instead of being infinite, the process takes the continuous 
spectrum of real values. For t > 1, the downwards jumps form a linear pure death 
process which "returns from 00" by a Poisson entrance law, again with the time 
transformed to run on the interval (1, 00). The deterministic time transformations 
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Figure 1. Simulated sample trajectory of the bi-Poisson process. 
The process follows the segments £j : y = 9{l — t)j — t/9, < t < 1, 
J = 0, 1, 2, . . . except for the upwards jumps, and then follows the 
half-lines y = 9{t—l)j—l/d, t > 1, j = . . . ,1,0 with the downwards 
jumps. 

are logarithmic and introduce a rather simple non-homogeneity into the birth rates 
and the death rates of the process. However, they force infinite number of jumps 
before and after t = 1. 

For a more formal description of (Xt), we set 







0<t<l 


(4) 


Xt = < 9Zi - i. 


t = i. 




[e{t-i)Zt-l 


t > 1, 



where random variables Zt are {0, 1,2,... }-valued for t ^ 1. We will construct the 
appropriate process {Zt)t>o in three steps: we first define {Zt)o<t<i as a pure birth 
process, then we extend it to i = 1 by passing to the limit, and finally we extend 
the process to <> 1 as a pure death process with a Zi-dependent Poisson entrance 
law. 

2.1. The pure birth phase. As {Zt)o<t<i we take the non-homogeneous linear 
pure birth process with immigration with the birth rate 

The properties of such a process are well known. For small enough \z\, jlUL Exercise 
5.1] gives the following generating function of the transition probabilities for the 
more general non- homogeneous linear pure birth process with the birth rate A„ (t) = 
v{t)+n\{t). 

^0 [l-z{l-pis,t))) 
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where p{s,t) = e~ ^^^'''^^ . In our setting, p{s,t) — jEj and with z = e"^-^~*^ we 
get 

where J^<s = a{Xr : r < s} = a{Zr : r < s}. Thus the conditional distribution 
£{Zt — Zs\Zs) is negative binomial with parameters r = Zg + l/O^, p = {l—t)/{l — s). 
(Tabled lists the parameterizations of the distributions we use in this note.) 



Name 


Parameters 


Distribution 


£;(e"^) 


Poisson 
Gamma 
Negative Binomial 
Binomial 


A > 
p > 0, fj > 

r>0, 0<j5<l 
n>0, 0<p<l 


e--'AVfc!,fc = 0,1,2,... 

^il^p^(l-p)^fc = o,l,... 

(^y(l-p)"-^ fc = 0,l,...,n. 


exp(A(e" - 1)) 
{l-au)-P 

(l-(l-p)e")'- 

(1 -p + pe")" 



Table 1. Laws that appear as transition probabilities. 



Differentiating |5l) at m = and using Q) we verify that {Xt, J-<t)o<t<i is a 
martingale. Setting s — 0, from © we get E{Xt) = and differentiating © again, 
after a calculation we get E{Xf) — t. 

In particular, Xi = lim4_,i Xt converges almost surely, as does Zi = limt^i(l — 
t)Zt. Taking the limit in Q we see that 

(6) E (e"^i I Z,) = (1 - u{l - s))-^^^-^/'' , 

thus C{Zi\Zs) is gamma with shape parameter p = + and scale parameter 
(T = 1 — s, see Tabled In particular, Zi is gamma with r — 1/9^, a — I, and the 
support of Xi is [— 1/0, oo). 

2.2. The pure death phase. We now extend (Zt) from 0<i<ltot>lby 
specifying {Zt)t>i as a pure death process with the death rate 

n 



t-V 



and the Zi-dependent entrance law C{Zt\Zi) which we take as the Poisson law with 
parameter A = — 1). Thus 

(7) E(e"^*|Zi) =cxpfzr''""^ 



1 



A well known property of the linear pure birth process is that for 1 < s < t 
the transition probabilities C{Zt\Zs) are binomial with parameters n = Zg, p = 
(s-l)/(t-l), so 

The Poisson distribution is indeed the entrance law: given 1 < s < t we have 

oo 

Pr(Zt = ^|Zl) = P^iZt = i\Zs = n) Pr(Z, = n\Zi). 
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Indeed, the right hand side is 

oo 

E 



^ (n - i)!(s - - 1)"-' 



We now verify that the two pieces of the process fit together into a well defined 
Markov process {Zt)t>o- For < s < f < i, by conditioning on Zi we get 

Pr(Z, = I E[zie-^^/^'-'^ 



j\it^iy{l-s)Zs+i/e-T{Zs + l/e^) 
T{i+j + l/9^) /t-l\ '+i/^' 



Thus is negative binomial with r = Zg + 1/9'^ and p = (i — f — s). In 

particular, Zt is negative binomial with r = 1/9'^ and p = 1 — 1/t. An elementary 
calculation shows that Xt defined by (^J has mean zero and variance t. 

A straightforward calculation leads now to the verification of the Chapman- 
Kolmogorov equations in the remaining two cases: 

(i) If < si < S2 < 1 < i and i,k>0 then 

oo 

Pr(Zt = k\Zs, ^''^^t = k\Zs, = j) Pr(Z,, = j\Z,, = i). 

Indeed, the right hand side is 

1 / ^-1 \ >^ T{j + fc + l/g^)(l ~ 52)'+'+^/'' {t - Iy-'{S2 - Si)'-' 

_ r(z + k + i/fl f t-1 y+^/'' (1 ~ / (t - i)(s2 - Si) ^^ 



A:!r(i + 1/6*2) Vl-si/ (i - S2)''+'''+i/'' V (l-si)(t-S2) 

^(^ + fc + l/02) /t_i y+i/e^ /l-si^'= 



fc!r(i + 1/6*2) \t-sij \t-si 
(ii) If < s < 1 < < i2 and i,k>0 then 

oo 

Pr(Zt, = k\Zs = i) = J2 ^ ^1^*1 ^ ^1^^ ^ 

n— 

Indeed, the right hand side is 

(l_g)fc(t,-i)Wi/e^ ^ r(n + i + l/g2) / (t2-ti)(l-.s) 

/fc!r(i + l/6i2)(ti-s)'+fe+i/e^(i2-l)fe^ (n-fc)! (^2 - l)(ii - s) 



r(A: + i + 1/6*2) 



i+l/e^ / , \ k 



kW{i + l/9-^) \t2-sj \t2 

Thus {Zt)t>o is a well defined Markov process which determines Markov process 
(Xt)t>o through the one-to-one transformation |0J. 
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3. Conditional moments 

We now verify that {Xt)t>n is a quadratic harness. 

Theorem 3.1. For 6 > 0, let (Zt) be the Markov process defined in previous 
Section. Let (Xt) be defined by Then (Xt) is the bi-Poisson process with 

parameters {0,9), i.e. it has covariance conditional moments and 0j with 
rj ~ 9 and q — 1. 

Proof. In SectionElwe aheady verified that E{Xt) = 0, EXf = t. Since C{Zt\Zs) 
is binomial for 1 < s < t, we have E(Zt| Zg) = j^Zg. Combining this with the 
aheady estabhshed martingale property for t < 1, we see that (Xt, J^<t)t>o is a 
martingale. From the martingale property we get (Q. 

To compute the conditional moments, we calculate explicitly the conditional 
distribution of C{Zt\Zs, Zu). These are routine calculations, so we just state the 
final answers, and omit most of the calculations of the corresponding moments. 

(i) If < s < t < u < 1 then C{Zt — Zs\Zs, Z^) is binomial with parameters 

n ^ Zu - Zs and p = [IZ'^^lu-l] - Therefore 

E{Zt\Ts „) ^Zs+ ~ ~ ""l {Zu - Z,) 

\ ti s,uj {\-t){u-sy ' 

_ {u-t){\-s) {\-u){t-s) 
{\-t){u-s) ' {\-t){u-s) 
Using (0J) we get 

?/ — t t — s 

^{Xt\ T,,u) = -tl9 + {X, + si 9) + {Xu + ulQ) , 

u — s u — s 

which gives Similarly, 
which gives 

Var {Xt\ T,,u) = ^^T'^^^I^^ (^(1 - s){Xu + u/9) - 9{l - u){Xs + s/9)) . 

{u — s)^ ^ ' 

A calculation gives ||2J). 

(ii) IfO<s<t<l<u then C{Zt — Zs\Zg, Z^) is negative binomial with 
parameters r = Zg + Zu + 1 /9'^ and p = . Therefore 

E{Zt\Ts,^.) = Zs + r{l -p)/p 
_ {1 - s){u - 1) {u-m-_sl ^ (u-l){t-s) 



{\-t){u-s) ' {\~t){u-s) " 9'^{\-t){u- sY 
which leads to |5J and 

r(l -p) 



N^x{Zt\T,.u) 



p2 



{u-l){l-s){t-s){u-t) ^^ 2^ 
= T, 7T27 \2 ^'^^ + Zu + 19 j, 

(1 — tY{u — sY 

which after a calculation leads to 
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(iii) IfO<s<l<f<u then C{Zt — Zu\Zs, Zu) is negative binomial with 

(t-l){u-s) 
it-s){u~l)- 



parameters r = + Z„ + 1/6'^ and p = . A calculation verifies 



© and 0. 

(iv) Ifl<s<i<u then C{Zt — Z^) is binomial with n = Zs — Z^ and 

P — (tZij(„Zsj ■ ^ calculation verifies and (j3Jl. 
The conditional moments for the remaining choices oi s < t < u follow by continuity. 

□ 

4. Additional Properties 

Proposition 4.1 (Poisson representation). Let (N^) and (M^) be two independent 
Poisson processes with intensity A > 0. If (Xt) is a bi-Poisson process with positive 
parameters (77, 9) then 

C (^(t (h{t)X_^ + i))^^^ , (^t (^Xs^ + i))^^J Xe/, = A - = £ {{N,^)t>o, {M^)t>o 
where 

Kt) = 
^ ' et 

Proof. Without loss of generality we assume rj = 0. By and the Markov prop- 
erty, it suffices to prove that 

(9) C = a) = £ {{N^)t>^) , 
and 

(10) C ^1 - a) £ ((Af,^)t>o) . 

Both equalities follow now from elementary calculations of finite dimensional dis- 
tributions using the conditional distributions identified in Sectional 

To prove iCOJl, take t„ < t„_i < ...ti so that 1 < h{ti) < h{t2) < ...h{tn). 
Then for fci > fc2 > ■ ■ • > fc„ denoting Yj — ^/i(tj) we have 

Pr(ri = h,Y2^k2,...,Yn = kn\Zi = A) 

= Pr(y„ = A;„|y„_i = fc„_i) . . .Pr(y2 - feln = A;i)Pr(yi = ki\Zi = A) 

A'^i exp(-Ati) ^^^^^ . Nfc„-i-fe„ .^fcl-fc2 



which proves pOfl . 

The proof of JHl is similar after using the generalized Bayes formula: for ti < t2 < 
...t„sothatO < l/h{ti) < \/h{t2) < ...l//i(t„) < 1. Then for fci < ^2 < • • • < fc„, 
denoting Yj = Zi/i^i^i.-^ we have 

Pr(yi = ki,Y2^k2,...,Y^^ kn\Zi = A) 

= "^""^'r-^t;^^^ Pr(^" = fc„|y„_i = fc„_i) . . . Pr(y2 = fc2|yi = fci) Pr(yi = fci), 

JZi (A) 

where fzi\Yn=k„ is the conditional density of Zi given y„ = fc„, and fzi is the den- 
sity of Zi , which are both gamma, see the last paragraph of Secti0n r2.ll Elementary 
calculations now prove □ 

Proposition 4.2 (Time-inversion). If (Xt) is a bi-Poisson process with parameters 
(0,0, q), then (iXi/f)t>o has the same distribution as {Xt)t>o- (Compare [HI, 



8 



WLODZIMIERZ BRYC AND JACEK WESOLOWSKI 



Proof. This follows from the fact that (tXi/t) satisfies (0, (0), and 0, and hence 
by 12 Theorem 4.5 and Proposition 4.13] is determined uniquely. For q = 1, the 
conclusion can also be derived directly from and the fact that Markov process 
(Zf)f>o has the same transition probabilities as (Z]^/f)i>o. D 

Proposition 4.3 (Distribution of upward jumps). For a bi-Poisson process (Xt) 
with parameters {9,9), where 9 > 0, define 

=snp{se [0,1) : Zs = i} , i^O,l,..., 

i.e. Ti is the time of the {i + l)-th jump of the process (Xt) from the line y — 
9(1 — t)i — t/9, < t < 1, for i = 0,1, . . ., see Fig. Q Then the joint density of the 
random vector (Tq, Fi, . . . , Ffc) is 

nn f ( F(^ + fc + i)(i-.,)A+^-^ 

(11) /(ro,r....,r.) (.0, .1, r (^) (1 - .oP(l - .r)^ . . . (1 - .._,)^ 

for < So < Si < • ■ • < Sfc < 1 (and otherwise). 



Proof. Let (Mt)t>o be a homogeneous pure birth process with birth rates A„ = 
n + 1/9^, n = 0,1, . . . . It is well known that the sojourn times tj of {Mt) in 
state j are exponential with parameter j + 1/9'^, j = 0, 1, ... . For < t < 1 



we have Zt = M_in(i_t), so -log(l - F^) = J2j=o'^3- Therefore u- = ln(l - 
Ffc_i) — ln(l — Ffc) (here, we set F_i — 0). Since the Jacobian of the transformation 



s 



j ^ ln(l - Sj_i) - ln(l - s.j),j = 0, . . . , fc is J(so, . . . , Sfc) = 11^=0(1 ~ ^o) ^ ^'^^ 
tq,ti, . . . ,Tk are independent, the joint density of (Fq, Fi, . . . , F^) is 

J(so, . . . , s,) n ({j + 1/9') exp (-{j + 1/9') In 

which simplifies to Hll|l . 

Alternatively, we can use the fact that 

R = P(so < Fo < si < Fi < S2 < . . . < Sk-i < Ffc_i < Sfc < F^) 

^ PiZ,,^0,Z,,^l,...,Z,,=k). 

The formula 

/(ro.ri,...,rfc)('So, si, ■ • • , sk) = (-1) t: — 7^ ^ — ^ — 

" dskdsk-i . . ■ dsidso 

yields Hll|) after a calculation. □ 

Proposition 4.4 (Distribution of downward jumps). Consider 

A, = inf{t > 1: Zt^i} , i = 0,l,..., 

i.e. Ai is the time of entrance of the process (Xt) onto the line y — 9(t — l)i — 1/9, 
t > 1, for i = 0,1, . . . . Then the joint density of the random vector (Aq, Ai . . . , Afc) 
is 

f f, , , , F(^+fc + l) it,-i)M''-' 

r(^)^r (io-i)^(ii-i)^..(t.-i-i)2 

for 1 < tk < . . . < ti < to (and otherwise). 
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Proof. From time-invertibility of the process, (Aq, Ai . . . , Afc) has the same distri- 
bution as (1/ro, l/Fi, . . . , l/Ffc). Thus 

/(Ao,Ai...,Afc)(ifc, ■ ■ • , ti,to) = -5-2 ;2 /(ro,ri,...,rfc)(l/*o, ■ ■ • , l/^fc), 

which simphfies to the expression above. □ 

Proposition 4.5 (Time to reach lower boundary). The time a bi-Poisson process 
(Xt) with parameters rj — 6 > 0, 9 = 1 reaches the horizontal line —1/9 on which 
it stays forever is finite but has infinite expectation. 

Proof. The distribution of Aq = mi{t > 1 : Xt = —1/9} is a special case of the 
distribution of jumps, but it is just as easy to derive it independently. Since Zt is 
negative binomial, for i > 1 we have Pr(Ao > i) = l-Pr(Zt = 0) = l-{l-l/ty/^\ 
From the inequalities (1 — x)p < 1 — px when 0<p<l,x>0 and (1 — x)p < 
(1 — x) when p>l,0<x<lwe get 



E(Ao) = / |l-(l-7) ]dt> min{l, 1/9^} 



dt 

— = 00. 

t 

□ 

Proposition 4.6 (Poisson limit). Forr]9 > let {xj:^'^"^) be the bi-Poisson process 
with parameters {ri,9), and let (Nt) be the Poisson process with parameter A = 1. 
As rj ^ the process {^xj:gf'') converges in D[0,oo) to the Poisson-type process 
{Nt~t)t>o. 

Proof. Calculating the conditional variances one can check that {^X^^f^) is a bi- 
Poisson process with parameters (r/^, 1) 

Consider now the bi-Poisson process (Xt) — {xj:^'') with parameters (e,e) for 
e = y/r]9. Then by the previous argument, process = {jXj:^l)f^Q has the same 

distribution as process {^xj:gf ^). Therefore it suffices to show that as £ ^ 0, the 
process converges in D[0, 00) to the Poisson-type process {Nt — t)t>o- 

We first verify the convergence of finite-dimensional distributions. For < ^ < 
1/e^, the appropriate version of Q is 

i41 = (i-te2)41-i, 

where Z^'^l is a (non- homogeneous) pure birth process on < t < l/e^. We 
will verify that the finite dimensional distributions of (Z^^l) converge to the finite 
dimensional distributions of (Nt). 

Fix arbitrary = to < ti < t2 < • • • < and ui, U2, ■ • ■ , < 0. It suffices to 
show that 



(12) lim E l^exp u^Z^f^^j j = [J cxp (-(i, - t,-i)(e^- - 1)) 

We rely on the following observation, which can be regarded as special case of 
Slutsky's theorem: if — > u and (Wi"\wt\---M"^) 

converges weakly to a 



10 



WLODZIMIERZ BRYC AND JACEK WESOLOWSKI 



random vector (VFi, W2, . . . , Wn) as a ^ and appropriate exponential moments 
exist, then 



(13) lim^E l^exp UjWj"'' + u^Wi'^^ | ) = E | exp | ^ UjWj + uWn 
From ||SJ) we have 




where 



and 



1 - eHn 



1 - sHn-i - e^{tn - t„-i)e"" 



This proves H12|l for n = 1, and shows that H12() holding for n — 1 implies (|12ll for 
n, ending the proof by induction. 

Therefore, the increments of (Z^'^l) are asymptotically independent Poisson- 
distributed with parameter A — t ~ s, and the finite-dimensional distributions of 
(Z^^l) converge to the corresponding distributions of {Nt). 

Tightness of in £'[0,00) follows from [7| Theorem 2.1 and Remark 4.2]. In- 
deed, their condition (i°°) holds as is a martingale and E{Y'^{t))'^ = t. Their 
condition (ii°°) holds as for fixed ai < s < 02 and all e small enough from Q we 
get 

Pr( at least two jumps in [s, s + S)) 
= Pr(X^^'' has at least two jumps in [se^, {s + S)s^)) 
< Pr(4%^, - Z^ll > 2) < 5'Ciaua2)Eil + e'ziff. 

Thus 

S^^ limsup sup Pr( at least two jumps in [s, s+6)) < 6C{ai, 02) — > as (5 ^ 0. 

e — *0 ai<s<a2 

Finally, their condition (iii) holds as for fixed 77 > and all a > small enough 
limsupPr(sup lY^s)] > rj) < limsupPr(Z^^) > 1) = ae"'' ^ as a ^ 0. 

£— *-0 s<a £-^0 

□ 

Proposition 4.7 (Brownian limit). Let {Xt)t>Q — {X^ )t>o be the bi-Poisson 
process with parameters ?] = 6. q = 1. Then as 6 the process (X^ )t>o 
converges in D[0,oo) to the standard Brownian motion {Bt)t>o- 

Proof. The convergence of finite dimensional distribution follows from uniqueness 
of the quadratic harnesses (see p.2 for the special case we need here) and the fact 

that the limiting process must satisfy (jJl by the uniform integrability of {(xj:^^)'^ : 
9 < 1}; in fact, E{X^) = 2t^ + {t + St^ + t^)9^. 

Tightness in D[0, 00) now follows from ^ Proposition 1.2], as (Xj^-* : <t < 00) 
is a martingale for each 9, and the limiting process is continuous. 
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Here we give a simple direct argument for the convergence of finite dimensional 
distributions. It is enough to prove that for all = to < ti < t2 < • • • < and all 
ui,U2, ■ ■ ■ ,Un close enough to zero 



(14) hm E l^cxp j j = ^ [TS*'^ - ^^-i) {j2 



We proceed by induction, suppressing 6 in to shorten the expressions. Without 
loss of generality we may assume ti < 1 (set ui — 0, if necessary). To verify H14|) 
for n — 1 and ti < 1 we use and with s — 0, which gives 

/ 1 f /I 



' = — r-a a -> exp -uth 

Suppose (23Il holds for = to < ii < ^2 < • • ■ < t;i-i < 1 and let t„ e (t„_i, 1). 
Then again using Q and (0) we get 

(15) E ^exp UjXt}j j = AeE |^exp + 
where 

, / l-t^ y/(e^(i-t„-.)) ^_«„(t„-t„_,) 

= : r. : ^„fl,.,,n-.,^ exp 



1 - t„-i - (tn - <„-i)e«-"(i-*") ; ^ ^ V ^(1 - tn-i) 

1 , N i/(e(i-t„-i)) / -, _ , 



' "',l-t„-i-(t„-t„_i)e«-(i-*")j '^Pr"l-^n-i 



Since 

Ao = 



1 . \ l/(9'(l-t„-i)) 



(1 - t„_i)ee«.(*.-*.-i) - (t„ - i„_^)een„(i-t„-i) ^ 

by Taylor expansion lim^^o — e2""(*"~*"-i). Similarly limg^o = u„, which 

using and induction assumption proves f(T^ . 

Now we know that (|14|l holds for all n > 1 and = <o < ^i < ^2 < • • • < ^n-i < 1, 

and we take t„ = 1. Using Q and 10 we see that p5|) holds with 

X i/(e(i-t„_i)) 

I e"" 

,1- (l-i„_i)u„0; 

and 

By (231 this proves (fTIjl . 

Now we know that 114|) holds for all = io < ti < i2 < • • • < tn-2 < tn-i = 1, 
71 > 1 and we take i„ > 1. By Q and Q we get ((T5|l with 

'ge«,.(t„-i) _ 

e(t„ - 1) 

and 

Ag = e("»~""^/^ ^ e^""(*"~^^ 
Again using ifT^ we get ifTHl . 



e"" = exp I — — I e"'' 
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Finally, we assume that (|14|l holds for all = to < ti < ^2 < • • • < tn-i, n > 1 
with tn~i > 1 and we take t„ > tn-i- 

Using Q and ® we see that (fT^ holds with 



tn — tn-l + (tn-1 — l)e' 
tn-l 



*u„t„-l 



l/(0(t„_i-l)) 



and 

By (231 this proves ((Tl|l . 



□ 
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